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Abstract 

We study by Monte Carlo simulations a model of knotted polymer ring adsorbing onto an 
impenetrable, attractive wall. The polymer is described by a self-avoiding polygon (SAP) on the 
cubic lattice. We find that the adsorption transition temperature, the crossover exponent <ft and the 
metric exponent v, are the same as in the model where the topology of the ring is unrestricted. By 
measuring the average length of the knotted portion of the ring we are able to show that adsorbed 
knots are localized. This knot localization transition is triggered by the adsorption transition but is 
accompanied by a less sharp variation of the exponent related to the degree of localization. Indeed, 
for a whole interval below the adsorption transition, one can not exclude a contiuous variation 
with temperature of this exponent. Deep into the adsorbed phase we are able to verify that knot 
localization is strong and well described in terms of the flat knot model. 



PACS numbers: 36.20.Ey, 64.60.Ak, 87.15.Aa, 02.10.Kn 



I. INTRODUCTION 



Like other forms of topological entanglement of polymeric chains, knots have relevant 
consequences for both physics and biology fl. It is known that they can be found in 
long closed macromolecules |2j, |3|, U, [5] , such as circular DNA [6| , and that they can affect 
important physical properties of them. For example, the migration velocity of circular 
DNA in gel electrophoresis depends on the knot type [7]. Among the various properties of 
knotted polymers, the determination of the length of the knots inside them, i.e. the length 
of the part of the chain which in some sense "contains" the entanglement responsible for the 
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ll|,ll2|. Indeed, 



overall knottedness, has attracted much attention in recent years 
the length of the knotted portion of the chain can be expected to play an important role in 
determining its physical properties. For example, the diffusion coefficient of knots tied into 
DNA by micromanipulation techniques should depend on the average size of the knots [3] . 
The action of topoisomerase on knotted DNA surely depends on how localized the knot is 
14j . The folding dynamics of a knotted protein 15] should also depend on the size of the 
knot. 



The determination, and even the definition, of the knot length, for real, 3D polymers, is 
however quite difficult. Thus, historically, the problem was first faced for simplified models 
in which one imagines that a knotted polymer ring is confined in two dimensions. In this 
way the conformations of the ring reduce to those of a polymer network in 2D. The network 
is made only of loops and, under simplifying assumptions, maintains a fixed topology. These 
objects are called flat knots HQ. 

One possible physical realization of flat knots, which allows to trace a link between them 
and "real" knots, is the following: think of a knotted polymer in 3D, and imagine it is fully 
adsorbed on an attractive, flat surface; the polymer will then become two dimensional, 
and will consist of loops, since the original 3D chain was topologically a circle. The loops 
are made of segments joined at vertices, which in general correspond to overlaps of the 
adsorbed 3D polymer on itself. Usually, in the flat knot model the information on the sign 
of crossings is not taken into account. One can study the length of the "knot" inside this 
object, at least under some simplifying assumptions: namely, the number of crossings must 
be kept constant and at the minimum value compatible with the corresponding topology 
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in 3D. In fact, in such a case the network will consist of a fixed number L of segments, 
while the knot length, £, can be unambiguously put equal to the total length of the L — 1 
smaller ones. Now, one basic question about the knot length is whether knots are localized 
or not. Knots are said to be localized if, sending the polymer length N to infinity, their 
(average) length (£) does not grow as fast as N. This means that, in the thermodynamic 
limit (N — > oo), the knot will behave as a point-like object with respect to the whole 
polymer. More precisely, the localization can be of two types: strong and weak. It is said 
to be strong when (£) grows slower than any power of N (e.g., as log(log(iV))), while it is 
weak if (£) ~ N l , with the exponent t strictly less than 1 (but larger than 0). When a knot 
is delocalized, {£) grows as fast as N: in this case the knot will always occupy an extended 
part of the entire chain. From the point of view of statistical mechanics, the determination 
of the localization behaviour of knots is a most interesting issue, since the exponent t is 
expected to be a universal (model independent) quantity. Thus, it is not surprising that 
research on knot length has focused on this aspect; this is true for flat knots as well as for 
3D knots, whose study inherited some terminology, and some ideas, of flat knot theory. 
For flat knots, Monte Carlo (MC) simulations, and theoretical calculations which employ 



the theory of polymer networks 



18], allow to make predictions on the value of t. It turns 



out that flat knots are strongly localized in the good solvent regime 



derealization transition, and become delocalized, below the 9 point 
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but undergo a 



The study of the localization behaviour of 3D knots is more recent, and has been per- 
formed employing original strategies and computer simulations for a consistent statistical 
definition of knot length [lo|, 12] . Indeed, the analytical treatment of the statistical 



mechanics of polymers constrained to have the topology of a nontrivial knot is very hard 



21] ; this is mainly due to the non-local character of the knottedness constraint, which makes 



impossible its description by a local hamiltonian, and thus prevents, e.g. , the use of standard 
field-theoretical techniques. The Monte Carlo results obtained in [lOL Il2| for (£) as a function 
of N have shown that prime knots in 3D are weakly localized, in the good solvent regime, 
with exponent t ~ 0.72. The weak localization of 3D knots and the value of t determined 



m 
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121 ] for SAPs on the cubic lattice have been subsequentially confirmed by simulations 



of off-lattice models 



Since flat knots should describe fully adsorbed knotted polymers, they are a useful 
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model per se, and not only for the indirect, qualitative insights they provide into 3D knots. 



Adsorbed polymers are in fact extensively studied 
is an important paradigm of polymer statistics 
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24| and the adsorption transition 



For polymer rings adsorbed on a plane it is known that knotting can occur, and it has 
even been pro ven, for specific models, that it occurs with probability 1 for infinitely long 



chains 



25|, 



26( . Thus, the behaviour of knots should be expected to play an important role 
in determining the physical properties of adsorbed polymers, as it does for swollen chains 
in 3D space. Flat knots represent only an extremely schematic model of adsorbed knotted 
polymers. In fact, a realistic model of an adsorbed polymer is given by a system consisting 
of a 3D polymer interacting with a short-range attractive, impenetrable, plane. It is known 
[27I that, in the case of unrestricted topology, this system exhibits a phase diagram with a 
desorbed phase at high temperature T and an adsorbed phase at low T; these two regimes 
are separated by a phase transition at a certain critical temperature T c . At T = the 
polymer becomes a fully adsorbed 2D object, except possibly for crossings (which are 
always present if the topology of the polymer is different from that of an unknot). Thus, 
the ground state of real adsorbed knotted polymers should be described by flat knots; 
but for any nonzero T, adsorbed polymers have excursions (i.e. connected bunches of 
desorbed monomers) which can be quite extended, even if they have finite length on average. 

So it is not clear if the true behaviour of real adsorbed knots for T > should be similar 
to that of polymer networks of fixed topology, or of their fully 3-dimensional counterparts, 



or, maybe, of something in between. Recent measurements 28], on samples of knotted DNA 



28| 



adsorbed on a substrate, indicate that knots are rather localized. However, the chains in 
seem not to behave as fully 2D objects. For example, the measured v exponent is between 
the 2D and the 3D value. Thus, a full explanation of the behaviour of these polymers may 
require to go beyond the simple flat knot model. The model of adsorbed polymers we are 
going to study may be useful in this respect. 

In this work we study the adsorption process of knotted ring polymers by means of Monte 
Carlo (MC) simulations. We focus on the adsorption transition and the low temperature 
(adsorbed) regime, for the simplest prime knots (3i, 4 1; 5i and 52). We check their ther- 
modynamic properties, in order to trace any significative difference between the behaviour 
of polymers with a fixed knot type and that of polymers with unrestricted topology. We 
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estimate the temperature dependence of the average knot length {£) and search for a possible 
transition between the T = oo regime, where knots are expected to be weakly localized, and 
the fully-adsorbed one (T = 0), which corresponds to flat knots and thus to an expected 
strong localization. 

The paper is organized as follows. In Section II we describe the model and the MC 
algorithms we use for the simulations. In section III we present the numerical results and 
discuss the knot localization properties in different regimes. We close this section with a 
discussion on the relation between knotted polymers in the strongly adsorbed regime and 
the model of flat knots. Section IV contains our conclusions. 



II. MODEL AND SIMULATION METHODS 



A flexible polymer ring of N monomers close to an impenetrable surface can be modelled 
by an A-step self-avoiding polygon (SAP) on the cubic lattice confined to the half-space 
z > and with at least one vertex anchored at the z = plane. To include a short-range 
attractive interaction between the surface and the polymer, an energy -1 is assigned to 
each vertex of the SAP having z = (visit). Denoting by v(u) the number of visits of a 
given configuration u the equilibrium properties of the model are described by the partition 
function 

z n (t) = ^y H/KsT (i) 

{«} 

where T is the absolute temperature and is Boltzmann's constant. If the sum in Eq. (pQ) 
extends to configurations u with all possile topologies, this model displays, in the thermody- 
namic limit, a second order phase transition from a desorbed (high T) phase to an adsorbed 
one (low T) 22|, |23j. In particular there exists T c > such that the limiting free energy 



T(T) = lim N' 1 log Z N (T) 

N^oo 



(2) 



is equal to — lo g K °, independent of T, for all T > T c and is strictly greater than — logA"° 



for all T <T r 



29(. The limiting value K° denotes the critical fugacity of standard, non 



interacting 3D SAPs 



30j. Let (v) be the average number of visits and 



p(T) = lim M 



(3) 
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the limiting fraction of visits. Then for all T > T c , p{T) = (desorbed phase) and for all 
T < T c , p(T) > (adsorbed phase). Right at the transition temperature T = T c « 3.497 



271 ] one expects 



p(T c ) ~ 



(4) 



where is the crossover exponent which is believed to be very close to 1/2 27j. Another 
way to detect the adsorption transition is by looking at the metric exponent v that controls 
the scaling of the radius of gyration through the power law 



(R g ) ~ N h 



(5) 



Indeed for T > T c one expectsz/ m 0.588, i.e the value for 3D SAPs 23j while for T < T c 
the value v = 3/4 of 2D SAPs [23J should be recovered. 



The adsor 
literature [22 



ution transition is one of the most succesfully studied transitions in the polymer 



231 ] . but if one restricts the sum of the partition function (CD) to SAPs with 



a given knot type very few results are available so far. One of these is the rigorous, strict, 
bound 

J*(T) < F(T) VT > (6) 



where JF®(T) is the limiting free energy (CQ) restricted to the set of unknotted SAPs 25l. 126] 
Inequality (EI) implies that for every finite va 



knotted goes to one as N goes to infinity 



25 



ue of T the probability that the polygon is 
261 ] . The only exception would be the set 



of SAPs lying completely on the plane z = (2D SAPs) i.e. the zero temperature limit 
for adsorbing rings. Of course, the knotting probability for finite N will depend on T and 
this has been investigated numerically 26(. No further studies have been performed so far 
on the effect of topological constraints on the adsorption transition but some reasonable 
assumptions can still be made. For example, one should expect that universal exponents 
such as the crossover exponent and the metric exponent v do not depend on topological 
constraints. On the other hand, the critical temperature T c , could change when switching 
from the unrestricted to a restricted topology ensemble, but there is no strong theoretical 
insight on what should happen. Also, we cannot tell whether T c should depend on the 
specific knot type or not. Given the difficulty in making any progress in this problem by 
analytical or rigorous means a natural way to gain insight is by Monte Carlo simulations. 
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For a fixed temperature T, SAPs with a fixed knot type are generated by using a Monte 
Carlo approach based on the BFACF algorithm, 31j. This is an algorithm which samples 
along a Markov chain in the configuration space of polygons of variable N and with fixed 
knot type. The statistical ensemble considered is thus grand canonical, with a fugacity K 
assigned to each polygon step. We adopt this algorithm because it preserves the topology 
and is irreducible within each set of configurations having the same knot type 32j. At a 
given T we used a multiple Markov chain (MMC) procedure 331 in which configurations 
are exchanged among ensembles having different step fugacities 34j. This is done in order 
to improve the efficiency of the sampling, especially at low T where the SAPs are strongly 
adsorbed on the plane. The temperatures we considered are T = 3.50 ~ T c and T = 1.25, 
2.00 and 2.75, all much less then T c . We also consider the value 1/T = (non interacting 
case) in order to compare it with the known situation of unweighted and geometrically 



unrestricted SAPs in 3D 



10 



121. 



Despite the use of the MMC sampling technique, the BFACF algorithm becomes quite 
inefficient in the strongly adsorbed phase and for high values of N. To improve the sampling 
in this regime we decided to use an hybrid scheme based on a combination of the BFACF 



with the pivot algorithm 



35j . Since any pivot move can change the knot type of the resulting 



SAP [36j, a check of its topology is needed before the move itself can be a ccep ted. This is 



37, 



done by calculating the Alexander polynomial A(z) in z = —1 and z = —2 

The BFACF algorithm does not preserve the value of N. Thus, in order to extract 
canonical averages at fixed N, we bin the data according to their N value. To collect 
enough statistics for a given iV we used bins of width 10. In this respect the symbol (..) 
indicates for us averages taken within a bin centered in N and with size 10. The knot types 
considered in the simulations are the prime knots 3i, 4i, 5i and 52- However, most of the 
results we present here refer to the trefoil knot (3i). 

For each sampled SAP the length i of the hosted knot r is measured by determining the 
shortest possible arc that contains the knot. The procedure works as follows [10. 1121] : given 
a knotted configuration, open arcs of different length are extracted by employing a recursive 
procedure. Each arc is then converted into a loop by joining its ends at infinity (i.e., at very 
far dinstance) with a suitable path. The presence of the original knot is finally checked by 
computing, on the resulting loop, the Alexander polynomial A(z) in z = — 1 and z = —2 
(see [l2 1 for details). In all the simulations considered we sample, for each value of K, over 
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20000 (independent) configurations. Since for fixed T a MMC scheme with 10 different K 
values is used, the total number of configurations considered in the statistics of a given T 
amounts to 2 x 10 5 . 



III. RESULTS 

A. Desorbed phase 

To check the validity of our approach we first compare the known situation of 3D swollen 
SAPs 0, Q with the one of SAPs confined in the upper half space by an impenetrable 
non attractive plane (non interacting case). Since the constraint z > should not play 
a significant role in knot localization, we expect strong similarities between the two cases. 
A first interesting issue concerns the the value of K c for the non interacting and confined 
problem compared to the one {K°) of the 3D case. For the whole class of SAPs (unrestricted 



topology) it is known that K c = K° 29] but no information is available about the relation 



between K c (t) and K°(t) i.e. for the corresponding quantities for SAPs with fixed knot t ype 



r. By using a MMC with 10 different K'b ranging from K = 0.2109 up to K = 0.2130 [39] 
we obtain a good evidence (wit 



that K c {t) = K° c {t) = K° c 



lin the confidence limit and for the prime knots considered) 
401 ]. We also confirm that at T = oo the metric exponent v 
coincides (within error bars) in the two cases and that is independent on the knot type 
(we estimate v = 0.59 ± 0.01 for all knots considered). We now turn out attention to the 
behavior of the average knot length (£) as a function of N. Previous studies have shown 
that for 3D swollen knotted SAPs 



(£) = AN 1 + o(N t ), (7) 



with t « 0.72 for the trefoil 



id, 



12J. Do we have the same behavior if the 3i SAPs are 
confined into one half space by an impenetrable (but still not attractive) plane ? This 
seems to be the case as witnessed by Fig. CE] where a log-log plot of (£) as a function of N 
is reported for the two situations. The two curves look indeed linear and parallel to each 
other, confirming a weak localization regime with an exponent t that is the same (within 
error bars) in the two cases. 
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FIG. 1: (Color online) Average knot length (£) as a function of N for SAPs in the bulk (bottom 
curve) and for knotted SAPs confined in the z > half space (top curve). In both cases the hosted 
knot is the trefoil knot (3i). (£) ~ N l holds in both cases with t = 0.73 ± 0.03. This value is 
consistent with the estimated for 3D SAPs in the bulk [101 . 1121 ] . 

B. Adsorbed phase 

When the attractive interaction between the SAP and the plane is switched on, the 
entropy is not any more the only ingredient in determining the equilibrium properties of the 
system. Instead, equilibrium is determined by the interplay between entropy and the energy 
gain in flattening the polymer on the plane. Moreover, when the topology of the ring is 
restricted to a fixed knot type, an additional entropic effect arises, because the entropy of 
rings with fixed topology and of rings with unrestricted topology are different. 

From a numerical point of view, simulations of SAPs in the adsorbed regime require more 
effort than those for the desorbed one. Let K c (T\t) be the critical fugacity for adsorbing 
SAPs at temperature T and with knot type r. For T >T C (desorbed regime), one can indeed 
assume, by analogy with the unrestricted case, KJT\t) = K c {r) = K®(r). On the other 
hand, for T < T c (adsorbed phase), K c {T\t) should decrease as the temperature decreases 
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4l| and one cannot rely any more on a known value of K c for simulations at a given T. 
Hence, for each value T < T c considered the value of K c (T\t) must be estimated first (by 
short MC runs) before collecting a significant amount of data for that temperature. 

Figure [2] shows the average knot size (£) as a function of iV for SAPs with a trefoil knot 
tied in. Different curves correspond to different temperatures ranging from T = 3.50, a value 
just above the adsorption transition for the unrestricted case, down to T = 1.25, a value deep 
into the adsorbed phase. The plot in the noninteracting case is also reported for comparison. 
One can notice that the N behavior of {£) for SAPs close to the adsorption point T = 3.5 
coincides with the one obtained for the non interacting case. Being the adsorption point 
the last point of the desorbed phase, this result shows that in this phase knots are weakly 
localized with exponent t that does not depend on T, and coincides with the one found for 
3D SAPs. Below the adsorption transition the situation changes significantly: the average 
knot length still follows a power law behavior (J7J) with N, but the exponent t decreases as 
T decreases. More interestingly, if we go deeper into the adsorbed phase (T = 1.5, 1.25) 
and look for sufficiently large N, (£) tends to a constant value. This is a signal of a strong 



localization regime (t = 0), reminiscent of the one found for flat knots [12]. 

Clearly if t = 0, Eq. ([7j) does not give any insight into the degree of localization of knots. 
More detailed information can be however obtained by analyzing the N behavior of the 



probability distribution function (PDF) of the knot length, P(£,N) [12J. In analogy with 



previous works on similar problems 
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431 ] one can assume the following scaling form: 



P(l } N) = r c -g(£/N 1J ), (8) 

where the scaling function g is expected to approach rapidly zero as soon as £ > N D , 
(D < 1). The quantity N D is a cutoff on the maximum value £ can assume. We expect 
D — 1, because there is no reason a priori to think that there exists some 'topological 
cutoff' which limits the size of the knot. This is confirmed by our measurements, which 
yield D ~ 0.9 4- 1 at every T. Assuming g is integrable when its argument is sufficiently 
large, one can deduce that, for < t < 1, c = 2 — t, while c > 2 always implies t = 0. In 
this respect the desorbed phase (where t ~ 0.75) is characterized by c = 1.25 1(3] , while at 
T = (fully adsorbed polymer) we could expect c = 2.69 i.e. the value found for flat knots 

m 

A common technique to analyze the scaling of the PDF (jSJ) goes as follows: for a trial 
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FIG. 2: (Color online) Average knot length {£) as a function of iV for SAPs with knot 3i. Different 
curves correspond to decreasing (from top to bottom) values of the temperature. The dashed curve 
corresponds to the non interacting case. The corresponding values of t can be deduced from the 
values of c reported in Table II. 

value of c and fixed N, one plots P(£,N)£ C versus E/N . Clearly, by varying N different 
curves are displayed but if the values of c and D are close to the correct ones, all these curves 
should collapse onto a single one described by g. Eventually, after several trials, "optimal" 
values of c and D can be estimated. Unfortunately, to have a good matching of the curves 
an extremely good statistics is required, and this would not be feasible in this context. 

We can instead perform an analysis based on the scaling behavior of the moments of the 
PDF in Eq. (1) 44j. This method relies on the following consideration: given the scaling 
behavior (jBJ) for the PDF, its g-th moment (q > 0) should obey the asymptotic law: 



(n 



£ q - c g(£/N D ) ~ ]\[Dq+D(l-c) = N t(q) 



(9) 



and the two parameters D and c can be deduced by fitting the estimated exponents t(q) 



11 



2.5 - 



1.5 - 




0.5 - 



0.02 
1/N 



0.04 



min 



FIG. 3: (Color online) Finite-size scaling analysis of the exponents t(q) for T = 1.25. For each value 
of q (q increases monotonically from bottom to top), t(q) is obtained by fitting the data of (£ q ) as a 
function of N with a power-law in the range between N m i n and N max = 1500, with N m i n <C N max . 
Here, different curves, corresponding to different values of q ranging from q = 1.00 (bottom curve) 
to q = 4.00 (top curve), are shown. To extract the asymptotic values we compute t{q) as a function 
of N m i n and extrapolate these values as 1/N m i n — > 0: This gives our best estimate of t(q). 



against the order q 



45| and performing a finite size scaling analysis (see Figure [3]). As 



shown in Figure H] for the T = 1.25 case, the plots of t(q) show deviations from linearity at 
relatively low q, due to finite N scaling correction effects. This is typical for this kind of 
analysis [44J. However, for a sufficiently wide range of q, a linear behavior can be identified 
whose intercept gives an estimate of c. Repeating the above procedure for the different 
temperatures considered one obtains the estimates plotted in Figure [5] and reported in Table 



For T > T c the value of c is roughly 1.25 independent on T. This is consistent with the 
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FIG. 4: (Color online) Analysis of the moments of the pdf P(i, N) for T = 1.25 and for the knot 
3i. The exponents t(q), calculated as described in the caption of the previous picture and in the 
main text, are plotted against q. For q > 2.5 a good linear behavior is obtained and a linear fit in 
that range of q gives the estimate of c. 



findings presented in the previous section. Indeed, in the desorbed regime (and up to the 
adsorption point included) the knot length exponent t is ~ 0.72. As the temperature is 
lowered the polymer goes more deeply into the adsorbed phase and the c exponent increases 
reaching, at T = 1.25, the value c = 2.55 ± 0.10. The curve in Figure [5] furnishes a good 
evidence that knots in the adsorbed regime become strongly localized. The value of c at 
T = 1.25 (2.55) agrees within the error bars with the value 2.69 found for flat knots, but it 
is possible that such value is reached precisely only in the T — > limit. 

One could take the intersection between the curve in Figure [5] and the line c = 2 as an 
estimate of a transition point T; oc between the weak and the strong localization regimes. This 
would suggest a localization transition occuring well below the adsorption transition. On the 
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FIG. 5: (Color online) Estimated values of the exponent c for SAPs at the adsorption transition 
(T = 3.5) and in the adsorbed phase (T < 3.5). The knot considered is the trefoil. For T > T c , we 
get c = 1.27 ± 0.03, independently on T (desorbed phase) This is consistent with c = 2 — t and the 
expected value (t rs 0.75) found in this regime. 

other hand, the curve in Figure [5] is an estimate of c(T) that relies on finite N simulations 
and it is hard to decide whether c(T) would show a sharp discontinuity as iV — > oo. Most 
intriguing would be the possibility of a range of temperatures in which the exponent c varies 
with T. 

Results for other prime knots are quite similar to those presented for the trefoil knot. 
Figure M shows for example the estimates of t(q) for 4i, 5i and 52 deep in the adsorbed phase 
(T = 1.25) compared with the one found for 3i. All the curves look quite similar and by 
performing a linear extrapolation we obtain the estimates of c given in table [III 

The estimated exponents are well compatible with each other. However, there seems to 
be a systematic difference between the behavior of 3i and 5i and that of the other knots, 
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T 


c 


3.50 


1.27 ±0.03 


2.75 


1.55 ±0.05 


2.00 


1.80 ±0.05 


1.50 


2.42 ±0.10 


1.25 


2.55 ±0.10 



TABLE I: Estimates of the exponent c for different values of T for trefoil knots. They have been 
obtained by the finite size scaling analysis of the moments of the knot length as explained in the 




FIG. 6: Estimates of the exponent t(q) as a function of q for T = 1.25 and for different prime 
knots. The error bars, that are not reported here for clarity, are of the same order as the ones 
reported for 3i in figure HI 
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knot type 


c 


3i 


2.55 ±0.10 


4i 


2.60 ±0.10 


5i 


2.55 ±0.10 


5 2 


2.56 ±0.16 



TABLE II: (color online) Estimates of the exponent c for different prime knots at T = 1.25. 

as can be seen from the plots in Figure [6j Note that 3i and 5i belong to the family of 
torus knots [37| and it could be that, for knots belonging to this family, (£ q ) displays an N 
behavior that is identical, not only in terms of the localization exponent t(q), but also in 
terms of the amplitude. The curve for the knot 4i is the most far apart and this is maybe 
due to the fact that 4i is achiral and belonging to a different family. 

C. Equilibrium behavior of the hosted knot. 

Having established that the adsorption transition drives knots from being weakly localized 
to be strongly localized, it is now interesting to understand whether the knot behaves as the 
rest of the chain. In fact one may wonder if the typical equilibrium configurations in the 
various phases are the ones in which the knotted part is expelled out of the plane, so that 
the knot is free to fluctuate in the bulk. This can be checked by comparing for example the 
average height (z) of the whole SAPs, (z), to the one restricted to its knotted part, (zknot)- 
Figure [7J shows the N dependence of (z), and (zknot), at two values of the temperature. At 
T = 3.5 (left panel) the two average heights are practically identical suggesting that, above 
the adsorption transition, the knotted part is indistinguishable from the hosting ring. This 
behavior seems to change at T = 1.25 (deep adsorbed phase), where (z) is systematically 
lower than (zknot) ■ This could indicate that in the adsorbed regime the knotted part tends to 
be, on average, further away from the adsorbing plane than the whole chain. Note however 
that, even in the strongly adsorbed phase, the knotted part must keep a minimal number 
of excursions, in order to connect the minimal number of crossings required by its topology. 
It turns out that this minimal number of excursions is sufficient to explain the differences 
shown in Fig. [7J (right panel). This can be seen as follows. By simulationing unknotted rings 
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FIG. 7: (Color online) Average height (z) of the SAP and of its knotted counterpart {(zknot)) as a 
function of N. The left panel refers to T = 3.5 (desorbed phase) while the right one to T = 1.25 
(strongly adsorbed phase). 

at T = 1.25 with N ~ 76 monomers (which is roughly the equilibrium length of the 3i knot 
at that temperature), we observe that the average number of monomers in the excursions is 
(b) « 32. On the other hand, the knotted portion of a knotted ring at the same temperature 
has (b) ~ 42. In both cases, almost all excursions have height z — 1. Hence the knot has, on 
average, ~ 10 more monomers in the bulk with respect to tis unknotted counterpart. This is 
in agreement with the discrepancy observed in Figure [7] (right panel) . This difference is then 
due only to the unavoidable crossings pertaining to the knot. In fact each crossing requires 
at least 3 excursions of lenght 3 each, resulting to a minimal excursion length 3 ■ 3 = 9 ~ 10 
for the knotted part. Similar considerations apply to different knot types. This suggests 
that the knotted parts of a polymer are essentially no more desorbed than the rest of the 
chain. 
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FIG. 8: Average number of visits (v) of the whole SAP (top curve) and its knotted part (bottom 
curve) as a function of N. The temperature considered is T = 3.50 ~ T c . The dashed line is still 
the one for the knotted one but now multiplied by {N/1)*. 

Given that the position of the knot within the chain and with respect to the adsorbing 
plane has nothing special if compared to any remaining part of the chain, it is reasonable to 
expect that the equilibrium properties themselves, once restricted to the knotted part, will 
display the same features as those of the whole chain. To confirm this picture we estimate the 
average number of visits as a function of N for the knotted part and compare it with the one 
of the whole polygon (Figured]). The temperature considered is the adsorption temperature 
T c where the scaling behavior (j3j) with « 1/2 is known to hold for all rings (unrestricted 
topology). By performing a simple power law fit we indeed get <fi = (0.53 ± 0.02) for the 
whole knotted SAP. Note that the estimate is for SAPs with fixed knot type, suggesting 
that the crossover exponent could be unaffected by topological constraints. To compute a 
crossover exponent restricted to the knotted part, the scaling H] must be replaced by 
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(Vknot) ~ (i)^, (10) 

where (vknot) indicates the number of visits of the knotted part of the ring. This gives 
(fiknot = (0.49 ± 0.04), compatible with = <p knot = 1/2. 

Similarly one can define the metric exponent of the knotted portion of the chain , Vknot 

as 

(Rknot) ~ W knot - (11) 

The estimates are reported in Table II I II and compared with those for the whole chain. 

The fknot determinations are always slightly higher than those of v and have slightly larger 
error bars. This is due to two effects: i varies on a smaller range than N, and the data for 
Rknot are noisier too. The measured values of Vknot are however comparable, within error 
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T 


V 




CO 


0.59 ±0.01 


0.61 ±0.03 


3.50 


0.60 ± 0.01 


0.64 ±0.03 


2.75 


0.75 ±0.01 


0.75 ±0.03 


2.00 


0.75 ± 0.02 


0.77 ±0.04 


1.50 


0.75 ± 0.02 


0.81 ±0.06 


1.25 


0.74 ± 0.02 


0.81 ±0.06 



TABLE III: Estimates of the exponent v for different values of T for trefoil knots. They have been 
obtained by a finite size scaling analysis of equations © and (fTTj) . 



bars, with the values of v for the whole SAP. At very low T, when knots become strongly 
localized the error bars on the estimates of Vknot are quite big. This is mainly due to the 
relatively small values of (£) (between 16 and 100) that do not allow a good asymptotic 
analysis of (TIT]) . 



D. Are adsorbed knots behaving as flat knots ? 

Our results show that the degree of localization of prime knots appear to be independent 
on the knot type (see table [Til)- This is consistent with the theory of flat knots which 
shows that, at leading order, all prime knots can be asymptotically described by the figure 
8 graph reported in Figure [TDl whose behavior in the good solvent regime determines their 



strong localization [17j]. However, adsorbed knots do not show the same degree of localization 
(measured by the exponent c presented above) as flat knots. In fact, they are always somehow 
less localized, i.e. they show values of c lower than that provided by flat knot theory. Note 
however that for T = 1.25 the c values for adsorbed and flat knots are comparable within 
2 standard deviations. On the other hand, since the Monte Carlo procedure deteriorates as 
T decreases it would not be useful to try to go deeper in the adsorbed phase. Indeed, there 
we where we would just estimate other values of c that would possibly match better with 
the flat knot value, but with larger error bars. 

We can then conclude that strongly adsorbed polymers should eventually behave as flat 
knots. If this happens at some T* with < T* < T c or just in the limit T — > it is however 
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FIG. 10: Sketch of the figure-8 graph. Here L is the total length of the ring while A is the length 



of the shorter loop, that can be identified with the knot length t. It can be shown [17[ that, when 
L >> A, the pdf for A scales as A _C8 , with c% = 2.69. 

impossible to establish within our numerical precision. 

Another distinction between adsorbed knotted rings and flat knots concerns the number 
of crossings, Cr, the chain makes with itself when it is regularly projected on the adsorbing 



plane 



461 ] . Indeed, in the flat knot model the number of crossings is minimal (3 for the 



3i, 4 for 4i, etc). This is crucial for the analytical treatment of the model and it would be 
interesting to see how close is the behaviour of the strongly adsorbed rings to this assumption. 

It turns out that while the average number of crossings, (Cr), of the whole chain 
increases (approximately linearly) with N, the same quantity restricted to the knotted 
portion is approximatively constant and close to the minimum value allowed by the 
topology. For example, (CV 3l ) 3.6. The independence of N is quite reasonable since deep 
into the adsorbed phase knots are strongly localized. Moreover, strongly adsorbed rings 
tend to minimize their energy by maximizing the number of visits and, in the limit T —>■ 0, 
we should expect all the unimportant crossings to disappear leaving a 2D ring with just the 
essential (topological) crossings defining the knot [47| |. 

At T ^ 0, however, fluctuations in the number of crossings are present and this could 
explain the deviation of the c exponent from the flat knot value for finite T values. In 
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this respect, it could be useful to study the behavior of flat knots when the restriction on 
the minimal crossing number is relaxed. Recently, Guitter and Orlandini [16] introduced 
a model of flat knots on a lattice where the ring is a 2D polygon with a number of 
crossings which can be changed and tuned with an appropriate crossing fugacity. By 
implementing our knot detection technique on this model we have estimated the average 
knot size for flat knots as a function of N for a wide range of values of the average 
crossing number. We find that flat knots are always strongly localized, independently 
of the average crossing number. This suggests that the crossing number is not a key 
feature to explain the value of c and that the flat knot model, even with a fluctuating 
number of crossings, is not a fair model for knotted rings in the adsorbed phase unless T — > 0. 



IV. CONCLUSIONS 

In this paper we study by Monte Carlo simulations the equilibrium properties of self 
avoiding polygons with fixed topology (knot) adsorbing onto an impenetrable wall. For un- 
restricted topology it is known that SAPs undergo an adsorption transition from a desorbed 
3D swollen phase to an adsorbed phase. We first show that for SAPs with fixed topology the 
adsorption transition is still present. Moreover, we give numerical evidence that the metric 
exponent v in the adsorbed and desorbed phases and the crossover exponent for SAPs with 
fixed topology, agree with the ones for the unrestricted topology case. Even non-universal 
quantities such as the critical adsorption temperature T c and the critical fugacity K° seem 
to be unaffected by the topological constraint. 

By using a novel algorithm that allows the identification of the knotted portion of the SAP 
[lo| we are able to focus on the equilibrium critical properties of this portion and compare 
them with those the whole SAP. We show that the knotted part behaves, in most respects, 
as any other connected subset of the ring with the same length. For example, we find that 
at T = T c the average energy of the knot scales as (£)^ where the value of <fi agrees with that 
of the whole ring. The metric exponent of the knot fknot-, which describes the scaling of the 
radius of gyration of the knotted part of the ring as a function of the knot length, is also 
consistent with the value of v of the whole chain at every temperature. Furthermore, the 
average displacement of the knot from the plane (height) is the same of that of the whole 
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ring, indicating that there is not a preferred height in space for the knotted portion. 

The main emphasis of our work is however on the localization behavior of the knotted 
portion of the ring as T varies. We find that, for T > T c , the average length (£) of the 
knotted portion grows as iV* where t ~ 0.72, consistent with the value found in [hj, [3] for 
knotted rings in the 3D bulk. This shows that knots in the desorbed phase and right at the 
adsorption transition are weakly localized, i.e. the presence of an attracting impenetrable 
plane does not change the localization properties of the knot. Below T c , the knot becomes 
more and more localized reaching a strong localization regime deep into the adsorbed phase. 
This crossover to more localized states is certainly triggered by the adsorption transition, but 
is quite smooth, as witnessed by the T dependence of the estimated localization exponent 
t. Thus, we can not exclude that below T c there is a continuous variation of the exponent 
t with T. A possible alternative scenario is the existence of a sharp localization transition 
at some TJ oc < T c , such that t has the bulk value for T > Ti oc , while t = (and c = 2.69)) 
below. To justify this latter scenario, one must assume that our data are affected by strong 
finite-size corrections, which cannot be numerically detected unless one performs simulations 
at much larger values of N. Note, however, that, even under this assumption, the possibility 
that T[ oc = T c should be discarded in account of the observation that t starts to decrease 
only below T c ; a crossover region is expected to be more symmetric around the transition 
temperature TJ oc . For sufficiently low values of T (deep in the adsorbed phase) the estimated 
c exponent agrees, within error bars, with the one found for flat knots. However, due to the 
large statistical uncertainty we cannot rule out the possibility that the flat knot regime is 
reached only in the limit T — ► 0. 

These results suggest that the relation between the flat knots model and adsorbed knots 
is nontrivial. The flat knot regime is reached only at low enough T (or possibly at T — > 
) where the number of crossings, arising from a regular projection on the adsorbing plane, 
is close to the minimal value dictated by the topology. On the other hand, in the whole 
adsorbed phase many excursions in the bulk are allowed and give rise to a number of crossings 
that exceed the minimal one and that can fluctuate widely as N increases. These fluctuations 
are certainly responsible for the deviations from flat-knot behaviour observed at high enough 
T. The flat knot model is certainly not an adequate representation of knotted rings in the 
whole adsorbed phase. 
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